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We investigate circumstances under which one can generalize Horndeski's most general scalar- 
tensor theory of gravity. Specifically we demonstrate that a nonlinear combination of purely kinetic 
gravity terms can give rise to an accelerating universe without the addition of extra propagating 
degrees of freedom on cosmological backgrounds, and exhibit self tuning to bring a large cosmological 
constant under control. This nonlinear approach leads to new properties that may be instructive 
for exploring the behaviors of gravity. 

I. INTRODUCTION 

Theories of gravitation on cosmic scales have become an area of intense interest, both as a possible explanation for 
the observed cosmic acceleration and as an exploration of consistent extensions of general relativity. Gencrically such 
extensions lead to additional degrees of freedom, e.g. scalar modes in scalar-tensor theories, with possible pitfalls of 
higher than second order derivative field equations that may lack a well posed initial value formulation, or of ghosts 
and other instabilities. 

Horndeski in 1974 wrote the most general scalar-tensor theory giving second order field equations in four dimen- 
sional spacetimes [1]. In an alternate view, Galilcon theories [2-4], shift symmetric scalar fields possessing nonlinear 
combinations of field derivatives, have recently been studied with interest as sound models capable of cosmic acceler- 
ation. Also recently, as an approach to solve the cosmological constant problem, a linear combination of four terms 
called the Fab Four has been identified [5-7] as the unique terms allowing self tuning vacua that can cancel a large 
bare A term. 

Here we draw on aspects of all three of these approaches to demonstrate that nonlinear combinations of terms 
involving shift symmetric scalar fields can possess interesting advantages and properties. This extension of the "most 
general" scalar-tensor theory retains second order field equations and avoids pathologies on symmetric spacetimes 
such as the usual Fricdmann-Lemaitre-Robertson- Walker (FLRW) and de Sitter cosmologies. While the nonlinear 
approach can be applied quite generally, we give a proof of principle using a simple example of purely kinetic couplings 
with noncanonical forms, extending the "purely kinetic gravity" of [8]. 

Our example employs a nonlinear combination of the standard kinetic term and derivative coupling to the Einstein 
tensor. Besides possessing second order field equations it does not add any further propagating degrees of freedom, 
and it can achieve lasting cosmic acceleration unlike the linear, canonical, purely kinetic gravity theory of [8], avoid 
at least some instabilities unlike the derivatively coupled Galilcon theory investigated by [9], and self tune away a 
cosmological constant like the Fab Four. 

Deeper implications exist beyond our simple proof of principle. We emphasize that the example given is intended 
purely as a proof of principle to inspire further investigation into the theoretical properties of general nonlinear 
combinations, and not as a fit to observations. 

In Sec. II we explain our nonlinear generalization procedure and the conditions under which no additional prop- 
agating degrees of freedom are generated. The equations of motion are solved in Sec. Ill on a FLRW background, 
giving the cosmic and field evolution complete with attractors, revealing two distinct ways of approaching a de Sitter 
asymptotic state. Section IV demonstrates the self tuning properties of the theory, erasing an initial cosmological 
constant. The perturbed equations in Sec. V yield the no-ghost and stability conditions and the evolution of the 
effective Newton's constant G c g. We discuss various implications of the results in Sec. VI. 



II. PROMOTION TO NONLINEAR FUNCTION 



The Einstcin-Hilbcrt action of general relativity is extremely simple, condensing all the gravitational influence into 
the Ricci scalar curvature R. To allow cosmic acceleration, however, one must add either a cosmological constant 
A or additional degrees of freedom such as a scalar field <fi, e.g. with a potential and canonical kinetic term X = 
(—l/2)gV- v <j) t j j (f) u where <f>^ = V M oV The cosmological constant, or the field potential, raises issues of fine tuning and 
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naturalness: why don't high energy radiative corrections affect the form and magnitude to something characteristic 
of the early universe? We therefore do not employ either (except in Sec. IV where we erase them). 

A canonical kinetic term cannot by itself give rise to acceleration but noncanonical (but still minimally coupled) 
kinetic terms can, called k-essence [10-12], or in the absence of any potential, purely kinetic k-essence (e.g. [13]). One 
can think of this as promoting the Lagrangian term linear in the canonical kinetic contribution to a function. This 
gcnerically gives an extra degree of freedom, in that the sound speed is no longer fixed to the speed of light. 

Alternately, one could promote the Ricci scalar term to a function, hence f(R) theories [14]. This again adds a 
degree of freedom and one can view these as coupled scalar-tensor theories. Similarly one can have theories involving 
the Gauss-Bonnet combination Gqb = R abca R a bcd — ^R ab R a b + R 2 of the Ricmann tensor, Ricci tensor, and Ricci 
scalar, either linearly in Gqb or promoted to a function [15, 16]. In generalized Galileon theories one can promote 
the coefficients of the standard Galileon terms to functions of the canonical kinetic term, for example [17, 18]. 

Taking these examples as motivation, suppose we take the Horndeski theory action, composed of the linear combi- 
nation of several terms, and instead promote them to nonlinear functions or combinations. In generality we cannot 
do this without resulting in higher than second order equations of motion or adding unconstrained degrees of free- 
dom. However, in specific circumstances we can. For example a theory involving a function of the Ricci scalar and 
Gauss-Bonnet term f(R, Gqb) can be sound [19]. In Fab Four terminology, this mixes George and Ringo (though we 
do not allow the field potentials). This is permissible because of particular symmetries within these terms. 

In purely kinetic gravity theories, similar symmetries impose a unique Lagrangian involving only the Einstein tensor 
coupled to the field derivatives [20] (which Fab Four term John basically replicates). Up to mass dimension 6, the 
action is just the linear combination of the canonical kinetic term and the Einstein coupled kinetic term [8] , effectively 
giving a disformal field theory. This could achieve transient cosmic acceleration but not an asymptotic de Sitter state, 
and was later shown to have ghosts [9]. Allowing for an arbitrary constant coefficient of the canonical kinetic term, 
one could achieve a de Sitter asymptote but the theory has an early time Laplace instability [9]. 

Merging these two approaches of nonlinear function promotion and purely kinetic terms of great simplicity, we 
examine as a specific example nonlinear functions of the canonical and the Einstein coupled kinetic terms. The 
combination of nonlincarity and noncanonical nature delivers new characteristics to the theory. Since this "hip-hop" 
kinetic evolution extends the Fab Four self tuning possibilities, among other properties, we call this new Lagrangian 
term Fab 5 Freddy. As the line "Fab 5 Freddy told me everybody's fly" from Blondie's Rapture [21] predicts, this 
term also enables cosmic acceleration and an asymptotic de Sitter behavior, indeed in multiple ways. 

The action we study in detail is 



where ff" is the Einstein tensor associated to the metric g^, S m represents the action for the matter fields, and M 
is a mass scale to keep cq dimensionless, where we normalize to M = Hq. When the function / is linear then this is 
the derivatively coupled Galileon (using only C 2 in [9]), generalizing the purely kinetic gravity model by allowing a 
free constant coefficient for the canonical term. 

To study the effects of the nonlinear promotion we consider two cases: 1) c\ = 0, so the canonical and Einstein cou- 
pled kinetic terms are directly coupled nonlinearly, and 2) c 2 = 0, so only the derivative coupling appears nonlincarly. 
This allows us to compare these two different theories with the same linear limit. 

We can rewrite the action in terms of a Lagrange multiplier field x, as 




^R + Cl X + f(c 2 X + ^G^^) + S, 



(1) 




R + Cl X + f(x) + (c 2 X + S^GT^ - X ) f- +S> 



(2) 



Varying the action (2) in terms of x we find 




(3) 



This has the solution 



X 



c 2 X + ^G^fr 



(4) 



except at particular points for which f xx = (and note that in the linear case x is moot). Subscripts x denote 
derivatives with respect to x- By re-inserting the solution Eq. (4) back into Eq. (2), we verify that we obtain the 
original action Eq. (1). 
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Introducing a Lagrange multiplier field \ helps understanding of the independent degrees of freedom. In particular, 
such a Lagrange multiplier can be coupled with other elements (such as the Ricci scalar, as in the f(R) theories, or 
the Einstein tensor, as in our case). Both the Ricci scalar and the Einstein tensor are functions of a second derivative 
of the metric. Therefore by integrating by parts, a time-derivative for the Lagrange multiplier may appear. In this 
case, such a Lagrange multiplier can in general acquire a kinetic term, and it may start propagating. This situation, 
as already said, is common to those theories which can be written in terms of a Lagrange multiplier coupled to a 
second-order operator, e.g. as in f(R) or R + /(Gqb)- The theory f(R, Gqb) introduces two Lagrange multipliers. 
This theory is quite interesting as it has been proven that only one of these two new scalar degrees of freedom will 
propagate on Fricdmann-Lcmaitre- Robertson- Walker backgrounds [19]. On the other hand, both these degrees of 
freedom do propagate on anisotropic backgrounds. Therefore whether or not these Lagrange multipliers propagate or 
not depends on the chosen theory. 

We will see that the theory at hand, Eq. (1), will not introduce on cosmological backgrounds any new degree of 
freedom. However, we will find that the high-A: limit (where k is the wavemode) dispersion relation of perturbations 
will be modified, leading to a scale-dependent speed of propagation, i.e. c 2 oc k 2 . This is indeed similar to what 
happens in the f(R, Gqb) case. On a formal level it will be interesting to study eventually our theory on anisotropic 
backgrounds to sec whether or not the Lagrange multiplier will start propagating and we will discuss this issue in a 
future project. 



III. EQUATIONS OF MOTION AND EVOLUTION 



We give the general covariant background equations of motions in Appendix A. Here we specialize to a homogeneous 
and isotropic spacctimc where the metric is FLRW. The theory then has the property that the equations of motion for 
the action remain second order. We include a barotropic fluid (i.e. matter and radiation) with energy density p and 
pressure P and assume spatial flatness for simplicity. The background equations of motion for the action of Eq. (2) 



arc then 



3Af p 2 ! H 2 = p + l -c^ 2 + \c 2 ] x 4> 2 + f xX -f + 9j^f x H 2 ^ 2 , (5) 



>p! " 172 fx?) H = -P-3 M 2 ,H 2 ^ f (6) 



\c^ 2 ~ 3 ^H 2 <j> 2 - X ~ 4 g w) f x + 2jf- 2 H<f> 2 f x , 



ci+c 2 f x + 6 1 ^H 2 f x ) 4> = -3 Cl Hj> -12-± f x j>HH - (f x + 3Hf x ) (c 2 + G^H 2 ) j> , (7) 



-^2 fx <PHH {fx + 3#/x) (ca + 6^ " 
X=(\c 2+ 3^H^ 2 , (8) 
p = -3H(p + P) , (9) 

where H = a/a is the Hubble expansion rate of the scale factor a. 

Linear perturbations about the background are important for calculating the growth of structure, which we consider 
in Sec. V, but also for analyzing the degrees of freedom. Details of the equations are given in Appendix B but here 
we note a key point. The coupled system of perturbed equations for the two metric potentials, the barotropic fluid 
density and velocity, the <p scalar field, and the Lagrange multiplier scalar field x does not possess any time derivatives 
8x- This indicates that x is merely an auxiliary field with no dynamics but rather an algebraic constraint, and is 
uniquely determined by the other fields. This arises because the Einstein tensor within fix) only depends on first 
derivatives and not second derivatives in the Robertson- Walker background, i.e. only H 2 appears. 

To obtain the solutions to the evolution of the expansion H and field <j>, we put the background equations in the 
form of an autonomous system of coupled equations, using the dimensionlcss parameters H = H / Hq and x = <fi' /M p i, 
where primes denote derivatives with respect to N = In a. Then 

x = \R ( 10 ) 

aa — \p 

H' = -l-^-x' (11) 
a a 
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where 

a = 2H- 6f x c G H 3 x 2 - 2c G f xx H 4 x 2 (c 2 Hx 2 + l2c G H 3 x 2 ) (12) 
p = -~2c G f xx il 4 x 2 (c 2 H 2 x + 6c G H 4 x) - Af x c G H 4 x (13) 

7 = ZH 2 + f x (|#V - X - 3c G ffV) +^+f+ ^H 2 x 2 (14) 

a = c 2 f x H 2 + 6/ x c G i? 4 + f xx H 2 x (c 2 + 6c G H 2 ) (c 2 H 2 x + &c G H 4 x) + Cl H 2 (15) 

A = f x c 2 Hx + l8f x c G H 3 x + f xx H 2 x (c 2 + 6c G H 2 ) (c 2 Hx 2 + 12c G H 3 x 2 ) + c x Rx (16) 

u = 3f x H 2 x (c 2 + 6c G H 2 ) + Zc x H 2 x (17) 

with Q r o the dimcnsionless radiation energy density today and 

X=|hV + 3cgHV. (18) 

To ensure the accuracy of our numerical solution we use as a check the constraint equation (5), written in the 
dimcnsionless parameters as 

ff2 = f^o + + 1 ^ /xx _ f + 9fxCGi}ix 2 + /x |fiV + |i/ 2 x 2 ) , (19) 

with f2 m o the dimensionless matter density today. The quantity in parentheses can be viewed as an effective dark 
energy density. An effective dark energy pressure can similarly be defined using Eq. (6), with the effective dark energy 
equation of state ratio = P^/p^. 



A. Early Time Evolution 



As in [9], one can identify the early and late time asymptotic solutions. At early times, during radiation or matter 
domination, when H 2 3> 1, one generally has \ ~ 3cgH 4 x 2 and fi^ < 1 (if one fine tunes the ci or c 2 terms to 
dominate instead then the energy density would decay as ~ a~ 6 and hence be uninteresting). In this case the 
solution becomes 

X ~ a 3[l+3J i . 6 +4e(l+u> 6 )]/[2(l+2e)] ^20) 

x ~ o -3(i-«*)/(x+a«0 5 (21) 

where Wf, is the barotropic equation of state (0 for matter domination, 1/3 for radiation domination), and e = xfxx/ fx- 
To go further we must adopt a specific form for the function /. Taking f(x) — Ax n , we have e = n — 1 and 

H „ -3n(l-« t )/(2»-l) (22) 

1 - n(l + w b ) 

W * = 2n-l ■ (23) 

While in the linear model (n = 1), the dark energy traces the matter during matter domination, this is not so in the 
nonlinear model. The case n = is a cosmological constant. Note that the dark energy is phantom (and has a ghost, 
we will later find) for < n < 1/2. This means that to avoid violation of early radiation/matter domination the 
field would have to be highly fine tuned, more so than a cosmological constant. For n ps 1/2, the evolution diverges 
(e = —1/2) and matter/radiation domination is violated. Note that this rules out functions that act like n = 1/2 
power laws at early times, such as a DBI type / = y/l + x — !• Therefore we concentrate on n > 1/2. 



B. Late Time Evolution to de Sitter State 



During its evolution, the model leads to cosmic acceleration near the present and an asymptotic de Sitter state. 
Interestingly, this can arise in two ways. For H ' = and also x' = as a fixed point one needs 7 = and either uj = 
or o = 0. Combining the expression for 7 with Eq. (19) leads to the condition 



H 2 x 2 [ Cl +f x (c 2 + 6c G H 2 )} =0 



(24) 
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This guarantees that u> = also. In the c\ = case, this is the same de Sitter point H 2 S = — C2/(6cq) as in [9] 
and exists irrespective of the functional form of f(x) ( as l° n g as / at the de Sitter point). Note that c 2 and 
cq must have opposite signs for this de Sitter point to be present. In the c 2 = case, there is a new de Sitter point 
H 2 dS = -ci/(6/ x c G ). 

Since for the c\ = de Sitter point we have \ — > 0, we should choose a function / such that /(0) = otherwise we 
are putting in a cosmological constant. But then no solution for this de Sitter point exists for the nonlinear power 
law model that simultaneously satisfies H — > const and p^ — > const ^ 0. Recall that f x ~ X™ -1 - 

However, there is yet another de Sitter solution that we can construct for our nonlinear model; this arises because 
of the evolution of x such that asymptotically x' ^ 0. This solution still has H 2 S = — c 2 /(6cg) but x ~ a 3 !™- 1 )/! 3 "- 1 ). 
Thus as H — > const, ir decays to for n < 1 while x diverges for n > 1. 

In Fig. 1 we exhibit H 2 and for the c\ = cases with n = 1.5 and n = 0.8. As noted, the n > 1 case grows 
quickly relative to the background components and so must start with a small (p^,/p m )i to preserve later matter 
domination. At a — 10~ 6 , say, this ratio must be less than 10 -6 ' 6 , but this is still not as severe as the cosmological 
constant fine tuning which requires 10 -16 . The n = 0.8 case can actually dominate over matter at a = 10~ 6 , but has 
rather drastic evolution at z rj 1 as it suddenly turns toward the de Sitter attractor. 




1+Z 



FIG. 1: The evolution of the effective dark energy density in the ci = case is plotted for power law functions / ~ \" with 
n = 0.8 (dotted red curve) and n = 1.5 (dashed blue). The expansion history (H 2 /Ho) 2 (solid black) is also shown (for n — 0.8 
though the n = 1.5 case is nearly identical on this scale). Note that for n < 1, typically must be set to large values initially 
so that it does not decay to too small values at late times; conversely, for n > 1 grows relative to matter and radiation and 
must be set to low values initially. 

For the c 2 = case, as mentioned we expect at early times no significant change to the dynamics since again the 
cq term will dominate over c\. At late times, since p§ contains terms with different powers of x there is no extra de 
Sitter solution (hence = const) with varying x. Thus the only de Sitter solution is ci + 6f x coH 2 = 0. Note that 
X freezes at a finite value and so / — > const. The energy density evolution looks quite similar to the c\ = case and 
so instead we show the evolution of the dark energy equation of state for the two cases in the first panel of Fig. 2. 
The spike in w^z) can be ameliorated by raising the initial field energy density, but as discussed above this would 
impinge on matter domination. 

Although both cases reach de Sitter attractors asymptotically, the manner in which they achieve this differs. For 
the nonlinearity only applying to the Einstein tensor coupled kinetic term (c 2 = case), the solution is the double 
fixed point H' = = x'. This holds as well for the full nonlinearity (applied to both kinetic terms, i.e. the c\ = 
case) when n > 1 (but this is not the stable attractor). In addition the full nonlinearity case also has a de Sitter 
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FIG. 2: Comparison of the cases c 2 = (solid black) and Ci = (dotted blue), both with n = 0.9 and [p t j,/p rn ](a = 10 -6 ) = 10 -3 , 
is plotted for the equation of state parameter (left panel) and the field evolution x — cj>'/M p \ (right panel) of the effective dark 
energy. The c 2 = case has a more extreme phantom feature in w(z) near the present, but both cases have the same early 
time tracking and late time de Sitter attractor (although for ci = one has x' 0). 



solution with H' = but x' ^ 0, i.e. 0^0. The field will decelerate, (j> < (accelerate, (j> > 0) for n < 1 (n > 1). 
The difference between the field evolutions for the two cases is shown in the second panel of Fig. 2, for n = 0.9. 

A summary of the de Sitter attractors is given in Table I, including results from Sec. V on the ghost and Laplace 
stability conditions of the perturbations. 



(ci,c 2 ) 


n 




ZdS 


XdS 


No-ghost Laplace 


(0,c 2 ) 


1/2 < n < 1 


-c 2 /(6cg) 








■/ ~ 


(0,ca) 


n > 1 


-c 2 /(6cg) 


CO 


00 


X X 


(ci,0) 


1/2 < n< 1 


-ci/(6c g / x ) 


const 


const 


/ ~ / 


(ci,0) 


n > 1 


-ci/(6c G / x ) 


const 


const 


X X 



TABLE I: Summary of de Sitter attractors is given for the two cases of the model, with a nonlinear function / ~ \ n - The two 
cases have different approaches to de Sitter, that merge in the common linear limit n = 1. 



IV. SELF TUNING 



The ability of the theory to reach a de Sitter asymptotic state without a cosmological constant is interesting, as is 
the overall expansion behavior of such a cosmological model, but more significant is the ability of Fab 5 Freddy to self 
tune, in the manner of John or Paul in the Fab Four [5]. This allows the scalar field <j> - even without a potential - to 
cancel an existing (large) cosmological constant. This even holds if the cosmological constant readjusts as it passes 
through phase transitions. Here we present a simplified analysis showing these key properties while neglecting matter 
or radiation components. 

The dynamical equations are identical to Eqs. (10)-(18) except for the replacement of fl r o a -4 by — 3J7a in the 7 
term, coming from (three times) the background pressure. (Note 7^ 0.7, the observed cosmological constant, but 
is instead the early universe, bare cosmological constant.) Two de Sitter points can be found, both of which are 
attractors. The first arises from the explicit cosmological constant, with H\ = 8itGpa/(3Hq) = Q\ and the second is 
the self tuning solution with i/f = — C2/(6cq), as we found in the absence of a cosmological constant. Note Hi 3> ^2- 

For the first solution, the scalar field contribution dies away as ~ x 2n ~ a -6n/(2n-i) ^ g0 ^ e p Ure cosmological 
constant is a fixed point of the dynamics. (Of course matter and radiation would also rcdshift away.) For the second 
solution the scalar field dynamically adjusts such that p<f, — > —pA- Note that unlike in the earlier sections p^ < 0. 
However the same approach to a de Sitter state occurs, with p^ dynamically canceling p\ and retaining a small 
positive residual energy density, evolving with x ~ a _3 ( n_1 )/( 3,l ~ 1 ) on approach to H\ —> —02/(600). 

These analytic behaviors are verified numerically in Figs. 3 and 4. We include a cosmological constant S1a = 10 8 
throughout the numerical calculation, and consider f(x) = — x" (adopting n = 1.5, cg = 1, ci = —5.6, ci = 0). We 
begin H away from both asymptotic solutions H = H± t 2, illustrating the behaviors for different initial conditions in 
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Figs. 3 and 4. 

In the left panel of Fig. 3 we observe the approach to the standard cosmological constant attractor with H 2 — > 
87rGpA/(3i?o); w hh the right panel showing the vanishing p^ —> 0. 

However, below a certain critical initial condition Hi (depending on the other parameters), we observe entirely 
different dynamical behavior. Now, H approaches the second asymptotic point H 2 = — c^HQcq). This occurs despite 
the large cosmological constant present in the model. We find that the absolute value of the <fi field energy density 
approaches the p^ ~ —p\ solution, canceling the vacuum energy in the field equations. Hence the model exhibits self 
tuning, for some range of initial conditions (e.g. H 2 (a = 1(U 6 ) 

To see how self tuning occurs, we must examine the equations of motion (here taking c\ = 0) 



< 10 6 H 2 for the parameters adopted in the figure). 



SM^H 2 = 



lc 2 fJ 2 + f xX -f + 9^f x H 2 t 



X 



;C2 



3^tf 2 
M 2 



(25) 
(26) 

(27) 



On-shell (that is, at the asymptotic de Sitter state), we have H 2 = — C2/(6cq) and H = 0, and hence the scalar field 
equation is trivially satisfied, carrying no information regarding the evolution of (f>. However, the scalar field equation 
contains an explicit a dependence, and the Hamiltonian density % in the Fricdmann equation retains (f> dependence 
on shell, both of which are conditions given in [6] for self tuning to occur. On approach to the de Sitter point, the 
scalar field continues to evolve while p^ and H approach constant values. 

If we choose initial conditions for Hi such that it is initially far from the attractor H 2 = — C2/(6cg), then one 
can use Eqs. (10) and (11) to calculate how the model approaches the de Sitter state. For the power law models 
fix) ~ x" '> we nn d that the dynamical behaviour of H is independent of n, and the evolution toward do Sitter has 
H — Has ~ a~ 3 , x ~ a 6 . Ultimately the evolution of H will depend on the presence of matter and radiation (which we 
have neglected here), and also the functional form of fix)- Whether a specific self tuning model can be constructed 
that gives rise to a viable cosmological evolution will be the subject of future work (see [7] for the Fab Four case). 





1e+06 



8itGpa/(3Ho) for high density 



FIG. 3: [Left panel] H evolves toward its standard cosmological constant attractor H 1 
initial conditions, here Hi = 10 4 y—cv/jGcG). [Right panel] Meanwhile the scalar field energy density decays away. 



Going further, we can verify that the self tuning also self adjusts if the vacuum energy undergoes a phase transition 
at some redshift. We numerically model such an energy density with a tanh function, and choose the pressure to solve 
the continuity equation P = — y — p. The evolution of the quantities H 2 , p^, and J7a are shown in Figs. 5 and 6, 
demonstrating that the two de Sitter solutions still hold and the self tuning mechanism remains effective. The explicit 
cosmological constant can be made effectively invisible in our model. 



V. LINEAR PERTURBATIONS 



Linear perturbations of the equations of motion are important for calculating the growth of structure and assessing 
the ghost-free and stability conditions of the theory. For subhorizon perturbations one adopts the quasistatic approx- 
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FIG. 4: As Fig. 3, but with initial conditions Hi = 10 3 a/— C2/(6cg). [Left panel] Now we observe the dynamics leads to the 
second attractor H 2 — > — C2/(6cg), despite the presence of a large cosmological constant. [Right panel] The scalar field energy 
density self tunes to cancel the pa contribution in the field equations. 




1e+06 
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FIG. 5: As Fig. 3 for Hi = 10 4 y— ^/(6cg), but with the large vacuum energy undergoing a phase transition. The standard 
attractor H 2 — > 87rGpA/(3i/o) applies, and the <f> field energy density asymptotically decays, subject to a mild jump at the 
phase transition. 



imation. We begin by using Eq. (B2) to write the equations of motion for the perturbation in the Newtonian gauge 
(j3 = in Appendix B) as follows: 



2A/ 2 V 2 $ = Pm 6 m + 9-^f xx H*^6 X 



CG 



'pi 



M 2 ' 



f x 2<^V 2 $ - lH<t>V*8<t> + f xxX Sx + ^fxx^Sx 



M 2 ' 



6 * = 2 w^ 2< * 



C2 



(28) 
(29) 



ci V 2 50 + c 2 / x V 2 ^ - c 2 <p(f xx 5 x )- + ^f x (2H + 3ff 2 )V 2 <50 - 6^H 2 ^(f xx S X J ~ c 2 f xx (4> + 3H<p)5 X (30) 
-G§2-fxxH 2 &x 6^f xx H(2H + ZH 2 )4>5 X - 4^/ x 0V 2 <f> - 4^/ x 0V 2 <f> + 4^/ x ff0(V 2 ^ - V 2 $) = 



0idj$ - + 5u(V 2 ^ - V 2 $ 



Af 2 



/ x {2(0 + H4>){ giJ V 2 54> - didjScfr) + 4> 2 [ 9ij 



V 2 $ + V : 



2H^ 2 (f xx S X ) ■ + f xx 4fl# + 2i?0 2 + 3tf 2 </> 2 ) S X ~ f xx X^X 



(31) 
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100 1000 10000 100000 1e+06 100000 1e+06 

1+z 1+z 

FIG. 6: As Fig. 5 with a large vacuum energy undergoing a phase transition, but with Hi = 10 3 ^/—c^JJEcq). Self tuning 
remains effective despite the transition and the second attractor H 2 — > — C2/(6cg) is approached. The (f> field dynamically 
adjusts energy density to cancel the vacuum energy, and the transition does not modify the Hubble parameter due to it already 
being on the attractor. 



There are some important differences between this case and the linear case where f xx = 0, i.e. no nonlinear 
mixing. Here, in the Einstein and (j> field equations, terms appear of the form fc 4< J>/ xx and k 2 &f xx , arising from 
S\- These will lead to scale dependence in the gravitational coupling strength G c g derived below. Recall that the 
standard Galileon case does not have scale dependent coupling on cosmic scales well above the Vainshtein scale (see 
Cos from [9]). 



A. Evolution of Gravity 



To investigate the modified Poisson equations defining the coupling of matter to the metric potentials, we can 
use the = i) perturbed Einstein equation to remove (f xx 5x)\ an d then substitute for the ^ i), (f>, and x 
equations. In the quasistatic limit appropriate for linear growth on subhorizon scales the (0, 0) perturbed Einstein 
equation becomes 



V 2 $ = ^a 2 G% ] Pm 5 m . 
The equivalent modified Poisson equations for the other metric potential combinations are 

V 2 </> = 47rG$Wm 



The gravitational couplings are 



G 



eff 



K 3 K 8 + 2« 2 K9 



Gn ki(k 3 k 8 + 2K 2 Kg) + K2(k 5 k s + 2k 2 «7 — ' 



(32) 

(33) 
(34) 

(35) 



G 



MO 

eff 



G 



Gn 

off 



1 eff 



K2 



Gn 



2k 8 



G 



(*) 

eff 



~ b eff 



K-2 



(36) 
(37) 



10 



where G { ^> = Gg } /G N and 

/ci = 1 - 124/ xx i/ 6 x 4 - c 2 c G f xx H 4 x 4 - c G f x H 2 x 2 (38) 

«2 = -2c G f x H 2 x (39) 

k 3 = ci + c 2 / x + 2c G / x (2-ffiT' + 3i? 2 ) (40) 

k 4 = -c 2 Hx - 6cG-ff 3 a; (41) 
k 5 = -Ac G f x H (Hx + H'x + Hx) - \2c 2 G j xx H b x 2 (3Hx' + IH'x + 3Hx) 

-6c 2 c G f xx H 3 x 2 (Hx + H'x + Hx) (42) 

k 6 = -c G f xx Hx [xk 2 + H (AHx' + 6H'x + 7Hx)] (43) 

*7 = 1 + c G f x H 2 x 2 - 2c 2 G f xx H 4 x 4 k 2 /a 2 (44) 

k s = -1 + c G f x H 2 x 2 (45) 

k 9 = 2c G f x H (Hx' + H'x + Hx) + 2c G f' x H 2 x (46) 



and all quantities are in dimensionless form, i.e. H = H/Hq, f = f /(M^Hq), and primes denote derivatives with 
respect to N = In a. 

In the de Sitter limit for the case c\ = 0, one finds G c s/Gn = However, although both x an( i x approach 

0, they do so such that f x \ — > 0, f x x 2 — > const, and f xx x 4 — > oo. Thus \kj\ 3> \ki\ — > +oo and G e g — > 0. That is, 
gravity appears to turn off at late times. This arises in this limit from the nonlinear structure of the theory, i.e. the 
presence of f xx and its power law behavior. 

The numerical solutions for the evolution G e g(z) are shown in Fig. 7. At high rcdshift the theory acts as general 
relativity, then deviations begin when (k/aH) 2 fl^ ~ 1. At this point, the K-? contribution to G e s will dominate due 

to the k 2 term. Since kj appears in the denominator of , this scale dependent effective Newton's constant will 
typically vanish at high rcdshift, during matter domination (the exact rcdshift will be scale dependent and will also 

be determined by the initial conditions for the scalar field energy density) . G^ on the other hand will not vanish at 
early times owing to its different K7 dependence. 

The gravitational coupling G^ entering matter growth behaves as GR until near the present, since large kj actu- 
ally cancels out from it. At low redshift it spikes and then vanishes. The gravitational coupling G^ + ^ entering light 

deflection is given by the mean [G^ + G^]/2 and so shows deviations at both high and low rcdshift. Wc empha- 
size that the current model is not proposed as observationally viable but rather to introduce interesting theoretical 
properties of nonlinear, noncanonical kinetic gravity. 



k=0.01 h Mpc 
k=0.1 h Mpc 
k=1 .0 h Mpc 1 



0.1 
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k=0.01 h Mpc 
k=0.1 h Mpc , 
k=1.0 h Mpc" 1 



100 1000 10000 100000 1e+06 
1+z 



FIG. 7: The gravitational couplings G c g in the nonlinear theory (left panel G^*\ right panel both for / = x°' 9 ) become 

scale dependent, as shown by the evolution for three different wave modes: k — 0.01 h/Mpc (solid black), k — 0.1/i/Mpc 
(dotted red) and k = 1.0/i/Mpc (dashed blue). Gravity vanishes at late times. 
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B. Ghost and Stability Conditions 



In order to find the ghost conditions, we need the action for the independent degrees of freedom. It is convenient 
for this task to evaluate Eq. (B2) in the flat gauge (i.e. $ = 0). Then we can see that the fields ip,/3 and Sx can be 
integrated out leaving only two scalars to propagate, i.e. 5<p (the new-gravity mode), and v (the matter mode). But 
there is a crucial subtlety: the quadratic term Sx 2 will generate a term proportional to k 4 S<j> 2 /a 4 . This means that 
this theory will modify the high k behaviour of the modes, and this will lead to possible cosmological signatures. This 
situation is similar to what happens for FLRW backgrounds in the J(R,,G G b) theories [19]. This fc 4 -dependent term 
vanishes when f xx = 0, that is when the action is linear in the combination c 2 A + (c G / 'A/ 2 )G? A "' '<f> ■ This behaviour 
is a typical signature of the presence of a massive mode {Sx), whose kinetic term vanishes, but not its mass. 

After removing the auxiliary field Sx, and this is possible only when f xx ^ 0, we can write down the action as 



A ab V a V b + ^ (diV a )(diV b ) - % (d 2 V a ) (d 2 V b ) - ^ (diV a )(diV b ) + Ce ab V a V b + M ab V a V b 



(47) 



where we have defined V\ — 8(f), V 2 — v. The matrices A,D,E,M, as well as the two coefficients B and C, are 
functions of the background. Here we have also defined e ab as the two dimensional antisymmetric matrix with 
(±2 = 1. Furthermore, the only non-zero matrix element of the matrix D corresponds to D\\. 
The no-ghost requirements are 

det[A] = [f xx {c G (c 2 + 6c G H 2 ) (l2c G H 2 - c 2 ) H 4 x 4 f x + H 2 x 2 (c 2 - 6c lC 2 G H 4 x 2 + 36c G H 4 + l2c 2 c G H 2 )} 

, \ "9 9 9 / -9 "9 9X 1 (1 + W)5 W (l — C G H 2 X 2 f x ) , N 

+c G (l8c G H 2 - c 2 ) H 2 x 2 f 2 + (c 2 + 6c G H 2 - Cl c G H 2 x 2 ) f x + cj — V 4A2 ^ > (48) 

A22 = [(4 (36c G H 2 + 5c 2 ) f x H 6 x 6 - c G H 4 x 4 (ac G H 2 x 2 + \%c G H 2 + 2c 2 )) f xx + 9c 2 G H 4 x 4 f 2 - 6c G H 2 x 2 f x + l] 

where in addition to the scalar field we have assumed the presence of a barotropic fluid with equation of state w and 



energy density p w = p w /( J ff^A/ 2 1 ), and 



A 2 = [{H 6 x 6 c 2 G (36c G H 2 + 5c 2 ) f x - c 2 G p w ii 4 x 4 - c G H 4 x 4 (c G aH 2 x 2 + l8c G H 2 + 2c 2 )} f xx 
+9c 2 G H 4 x 4 f 2 ~ 6c G H 2 x 2 f x + l]w- c 2 G p w H 4 x 4 f xx . (50) 

During the radiation era, A 2 « w = 1/3 and the square brackets in A22 resolve to 1, so indeed A22 > 0. In the 
matter era where w = 0, then A22 has the same sign as —fxx an d so we require n < 1 in the power law model / ~ x n ■ 
For det A, since H 2 ^> 1 then the cq terms will dominate in the early universe over the other scalar field terms in the 
absence of fine tuning them to be small. This results in the condition (2n — 1)/(1 — n) > 0, satisfied for 1/2 < n < 1. 
Thus such theories are free of ghosts. 

Checking the speed of propagation of the field, with Laplace stability given by nonnegative sound speed squared, 
c 2 > 0, is somewhat more involved. In the high-A: limit, we find that the dispersion relations are given by 



2 = B 2 + A 2 2Dn k 4 _ 164/ 2 ffV7 xx k 4 



2 D11E22 k 2 k 2 

1 B z + A22D11 a 2 - a z 



where A is defined as 



A = Sxk [cg(c 2 + 6 c G H 2 )(c 2 - 12 c G H 2 )H 4 x 4 f x + 6 Cl c 2 G H 6 x 4 - (c 2 + 6 c G H 2 ) 2 H 2 x 2 ] 

f 2 + [c lCG H 2 x 2 - (c 2 + 6c G H 2 )]f x 

The speeds of propagation arc then found as the group velocity c = adu/dk, or 



c G (l8c G H 2 - c 2 )H 2 x 2 f 2 + [c lCG H 2 x 2 - (c 2 + 6 c G H 2 )]f x - c x . (53) 



64c 4 f 2 H 8 x 6 fxx k 2 9 
4= n \ ^ 2 f A ->0, c 2 pl = w>0. 54 
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One of the two solutions is trivial as it corresponds to the speed of the perfect fluid, but the other one sets a stability 
condition, and states that the speed of propagation will be scale dependent. 

Because of the fc 4 terms in Eq. (47), the dispersion relation and hence sound speed will be wavenumber dependent. 
From Eq. (54), we see that in the high k limit the leading order contribution to c 2 goes like 

(55) 



where CIa = SirGp^/ (3H 2 ). During matter domination, we typically find k/(aH) ^> 1 for sub-horizon modes relevant 
to linear perturbation theory, and f2^ ~ O (10 Q m J (this is largely dependent upon the initial conditions imposed, 
however this is a conservative upper bound on how large can be during matter domination). Hence we expect the 
fc 4 term in Eq. (47) to be the dominant contribution to between z ~ (1, 1000). 

If this held for radiation domination, then c 2 ~ (k/aH) 2 ^ (1 — n)/(2n — 1) and so the theory would be Laplace 
stable for 1/2 < n < 1. However, in the linear case n = 1, the terms proportional to f xx vanish identically and the 
leading order k 2 contribution vanishes, leaving a scale independent sound speed. As found in [9], the linear theory 
composed of a standard kinetic term and a kinetic term coupled to the Einstein tensor (i.e. the purely kinetic gravity 
theory of [8] with the canonical kinetic term generalized to have an arbitrary constant coefficient) is Laplace unstable 
in the radiation era. Thus, the nonlinearity of the current model can avoid that instability. 

However, at early enough times during radiation domination fl^ drops so low that the (k/aH) 2 ft{ term becomes 
subdominant for the modes k relevant to linear perturbation theory. To calculate the leading order behaviour of 
in the very early Universe, it is more instructive to consider the issue from a different angle. 

To preserve the CMB acoustic peak structure, and also to obtain an expansion history consistent with observations, 
we want initial conditions during radiation domination such that the effect of the scalar field 4> on the background 
expansion and the metric perturbations is negligible. In this case we can assume that fi^, <C 1 and the metric 
potentials are sourced by density perturbations only, hence we can treat the scalar field perturbations 8<j) as evolving 
on an otherwise standard cosmological background. Under this assumption, one can analytically calculate the no- 
ghost and Laplace conditions. For the linear model f(x) ~ Xj it wa s found in [9] that the scalar field perturbations 
possessed a sound speed c| that was negative during radiation domination; A = —1/3. This would lead to exponential 
growth of the scalar field perturbations, destroying the standard cosmological picture. One can perform a similar 
analysis for the more general f(x) case. Wc find the following no-ghost and Laplace stability conditions 

3c G ^ 2 (f x + 2 x f xx ) > (56) 

and 

2H + 3H 2 f x 



3H 2 f x + 2 X f- 



> (57) 



xx 

where we assume that at early times the cq contribution to the <fi energy density is much larger than the standard 
canonical C2 term (which is valid barring an extreme fine tuning of cg). For the power law models, during radiation 
domination these conditions correspond to 

c G /x < (58) 
2n - 1 < , (59) 

which would violate the positivity of p$. This could potentially cause problems in the transition to the matter 
dominated era, where p$ > is required to ensure that the no-ghost condition is satisfied. We conclude that the 
power law models cannot simultaneously satisfy the no-ghost and Laplace stability requirements at all times while 
also having p^ > during radiation domination. This does not preclude the possibility that a non-power law model 
might be constructed that can. 

The presence of a scale-dependent speed of propagation is a feature of this model and it has physical implications, 
especially at late times when f2^ — > 1. The reason for the presence of such a term may be due to the large symmetries 
of the FLRW manifolds, similarly to what happens in the context of f(R, Gqb) theory (see e.g. [19]). In that case it 
was shown that the kinetic term of one of the scalar perturbation modes was vanishing on general FLRW manifolds, so 
that it could be integrated out from the Lagrangian, giving rise in this way to a scale dependent speed of propagation 
c 2 cx k 2 /a 2 . Furthermore, also in that theory, at late times as the fc 2 -regime starts dominating, gravity for high k's 
tends to become weaker and weaker, i.e. G e ff/G/v — > 0. It would be interesting to study, along the same lines of the 
/(-R,Ggb) theories, whether anisotropic models (such as Bianchi-I type manifolds) possess more propagating degrees 
of freedom than FLRW. This will be investigated in a future project. Nonetheless, if indeed this scenario does happen, 
then this theory would behave similarly also to massive gravity, as it was shown that for that theory the kinetic terms 
of three perturbation modes vanish on FLRW due to the high degree of symmetries of the background [22]. 
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VI. CONCLUSIONS 



Gravitation is a fundamental force that we have just begun to explore cosmologically. One of the great advances 
made in gravity research in the past few years is the realization that symmetry principles both strongly restrict 
the theory and open up new avenues and effects. Galileon gravity and massive gravity both use shift symmetric 
fields and their couplings to functions of the metric to enable new properties, including cosmic acceleration without 
a cosmological constant or field potential. An action allowed by the symmetries and well behaved in initial value 
formulation, specifically one leading to second order equations of motion, is of particular interest. If moreover the 
field exhibits self tuning, allowing it to overcome a high energy cosmological constant, the theory is well worth 
examining. 

We show that by promoting a purely kinetic gravity term to a nonlinear function, possibly mixed with a noncanonical 
kinetic term for the field, fascinating properties can ensue. In addition to second order equations of motion and self 
tuning, the theory does not incur extra propagating degrees of freedom on a cosmological background. Similar effects 
of symmetric backgrounds are seen in massive gravity. For example, in massive gravity it was found [22] that in 
isotropic spacetimes the shift symmetric (Stiickelberg) fields have vanishing mixing between the graviton and the 
scalar mode, and furthermore their kinetic terms vanish. 

The new term discussed here, "Fab 5 Freddy," can self accelerate and is merely the harbinger of a whole class of 
such nonlinear promotions or combinations of terms. 

The background evolution of the expansion and field lead to early time tracker behavior and late time de Sitter 
attractors. Solving the linear perturbation equations we see that simple power law functions can be free of ghosts. The 
gravitational coupling and dispersion relation of perturbations become scale dependent, possibly leading to an early 
time instability and a late time vanishing of gravity. The specific models studied may not be observationally viable 
but the characteristics arising from the nonlinear, noncanonical action open new aspects of gravity. Most intriguing 
is the self tuning property that can cancel a bare cosmological constant dynamically, even through phase transitions. 
The evolution of the field basically makes A invisible. 

That the most general scalar-tensor theory giving second order field equations in 4D could be further generalized, at 
least on cosmological backgrounds, is highly interesting. The specific term considered here, a nonlinear promotion of 
the field kinetics coupled to the Einstein tensor (the unique, low mass dimension shift symmetric combination giving 
second order field equations), is merely a proof of principle, while theoretically instructive. Ways to extend this class 
of theory more generally, to different nonlinear functions and combinations, are straightforward and may preserve the 
most interesting and desirable characteristics while leading to more viable predictions experimentally. 
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Appendix A: Covariant Equations of Motion 



The scalar field equation is given by 



Cl D0 + c 2 V Q [f x V a 4>] - 2c G G^V Q [f x V 4>] = . 



(Al) 



The x field is given by 



X = 



■■ p V a <i>Vp<j> . 



(A2) 



The Einstein equations are given by 



G 




(A3) 
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where 



-cg/x 



+2c G 



1 



+V Q / x V Q ^V AI V„0] + s^„/ - s M „/ x x + ci 
and parentheses in a subscript denote symmetrization of the indices 



+ C 2/x 



(A4) 



Appendix B: Perturbation Equations in Detail 
1. The scalar perturbations 

Let us write down the perturbed metric in the following form 

ds 2 = -(1 + 2^)dt 2 + 2d i j3dtdx i + a 2 (1 - 2$)cte 2 . 



(Bl) 



Expanding the scalar field as <j) = <p(t) + 5(j>, and considering a barotropic perfect fluid with equation of state P = wp 
(for an action approach of perfect fluids see e.g. [23]), then we find that in Fourier space, the action at second order 
in the perturbation fields can be written as 

= J dtd 3 x a 3 j- (Witp + W 2 S(t> - W z k> - W±Scj) - p (1 + w) V + W 5 S X ) ^jr + \ (^~~~ _ w *j ^ 2 



w 7 



1 



1 



(1 + w) (V 



V -3wHV 
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w 



--W- 



w 9 



d 2 Scj) 



<9 2 $ 



1 



,-,2 1 



W w — + WuSx \if> + r W 12 5<j) + t: W 13 $ 



16 



(<9$) 2 lp(l + w)F 2 lp(l + w)k 2 V 2 1 



2 fxx 5 X 



WnS(f> - W l8 $ + W 19 ^ r )6 x - l-W 2 o-^--9pwH(l + w)V + 3 (1 + w) P V + W 2 i8cf> J $ L (B2) 

where the matter field is the scalar component of <5T°; = — p(l + u/JcJjV, so that the matter density contrast 
S m = Sp/p can be written as wS m /(l + w) = V — 'SwHV — ip. 

Notice we still have one gauge degree of freedom to choose. For example, we can consistently set fj = (Newtonian 
gauge), or 8<t> = (uniform field gauge), or $ = (flat gauge). 
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The coefficients of the previous action are the following: 
W g 

m - 

M 2 



W 2 
W 3 
Wi 
W 5 
W 6 

w 7 

W 8 
W n 
W 12 

Wis 
W u 
W 15 
W 16 
W 17 
Wis 

Wl 9 
W 20 
W21 



2HM 2 1 -6H^f x 4> 2 , 
4>c 2 f x + <t>c x + 6 ^Y^2-H 2 f x , 



= -citf + 6 M 2 ,H 2 - 36 ^f x H 



c 2 f x 4> 2 



l%^H 2 f x + ^c 1+ ^c 2 f x 
6 H A/ 2 - 18 H^fJ 2 , 



2 c 2 + 12f x 



f xx cp c 2 f i*J xx <pi m2 



2CG ^ 



-6M p V 

CG 

A/ 2 ' 



M 2 



c 2 f x + 4j^f x H + c 1+ 6c G f x H 2 , 



A/ 2 ^ 



il/ 2 



M 2 



2 A 



6 2 ^ 
Af 2 ' 



4 l|/x-4^ 2 (6 CG i? 2 + c 2 )/ xx 



Af 2 ' 



Af 2 



{Q^fxH 2 + c 1+ c 2 f, 



4 CG/y/0- 24^ / xx 



The equations of motion for the perturbations in any gauge can be derived by using standard 



(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 
(B15) 
(B16) 
(B17) 
(B18) 
(B19) 

b 3 HH , (B20) 
(B21) 

variational calculus. 



2. Tensor perturbations 



By introducing the two polarizations of transverse and traceless perturbations for the metric, we can write down 
the action expanded at second order as 



S GW = E / 



'pi ~M2 JxV ) x 



{dh x f 



pi 



^- f 6 2 ) ( — 



(B22) 



This gives the no-ghost condition 1 — cqJ x H 2 x 2 > 0, and speed of propagation equal to 



C GW — 



1 + c G f x H 2 x 2 
1 - CGf x B 2 X 2 



(B23) 
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Note that due to the coupling to the Einstein tensor this is not equal to the speed of light. A stable evolution for the 
background requires that Cq W > 0. 
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